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Abstract 



Q^' We present the next-to-next-to-leading order relation between the moments 

. of the T system spectral density and the inclusive S-meson semileptonic width. 

- ^ , The perturbative series for the width as an explicit function of the moments is 

' well convergent in three consequent orders in the strong coupling constant that 

provides solid and accurate theoretical estimate. As a result, the uncertainty 
of the value of \ Vcb\ Cabibbo-Kobayashi-Maskawa matrix element is reduced. 



The inclusive 5-meson semileptonic width rsi(6 clvi) is rather a clean place to 
obtain the value of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |Vcf,| 
(see refs.[|l], |], |] as a review). The main uncertainty of this estimate is related 
to the strong dependence of the result on the 6-quark pole mass which has to 
be known with extremely high accuracy to get a precise estimate of \Vch\- On the 
other hand the high moment sum rules for the system of T resonances are very 
sensitive to nib and can be used for the precise determination of this parameter after 
the proper treatment of Coulomb effects 0, §]. In the present paper we construct 
the direct relation between the moments of the T system spectral density and the 
inclusive i?-meson semileptonic width up to the next-to-next-to-leading (NNL) order 
of perturbative expansion in a^. In this way we avoid the strong dependence of Fgi 
on mb and essentially reduce the theoretical uncertainty in |Vc6|. It is known that 
if the moments for the T system and the inclusive semileptonic width are expressed 
through the pole mass mf, the perturbative expansions in for both quantities seem 
to be divergent. The heurictic criterium of convergence of perturbative series, i.e. the 
requirement that next correction is much smaller than the previous one, is violated. 
However, the explicit dependence of the inclusive semileptonic width on the pole mass 
TTib can be removed by reexpressing it through the moments. As a result one arrives 
at the perturbative expansion for the quantity Vcb which converges well up to NNL 
order in a^. 

The theoretical expression for the inclusive semileptonic width of B-meson reads 



where F{x) = 1 — Sx — 12x^lnx + — x^. The perturbative coefficient Cr(«s) is 
known up to the second order in 
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From formula (P one sees that the semileptonic width depends rather strongly on 
mi,, the b-quark mass, which is the pole mass. This means that if rrib is taken from 
some other experiment it should be determined with great accuracy for reasonable 
determination of \Vcb\- Being analysed independently, the perturbative series in 
for physical observables expressed in terms of the pole mass m?, seem not to enjoy a 
fast explicit convergence that can lead to large uncertainty of the numerical value of 
rrih extracted from different experiments. However rrif, is not an observable and has 
no immediate physical meaning. Therefore it can be safely removed from relations 
between physical observables. Here we use this strategy for direct determination 
of the mixing angle \Vcb\ reducing theoretical uncertainties to a great extent. Our 
analysis consists in direct relating the factor ml in eq. ([l|) to the fifth moment of the 
T sum rules 



The moments are defined as normalized derivatives of the 6-quark vector current 
polarization function H(s) 



and R{s) = 127rImH(s+2e). The (dimensionful) experimental moments are gen- 
erated by the function Rb{s) which is the normalized cross section Rb{s) = cr(e+e^ 
hadronSft5)/cr(e+e~ — > fJ^^fi^) 




(2) 




where 





(3) 




(4) 
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The theoretical moment Ai^^ is a dimensionless quantity which depends on 



Here Qb = —1/3 is the 6-quark electric charge. 

only logarithmically as In (mb//i) (in a finite order in as)- Thus, the substitution of 
relation (|]) to eq. (P substantially reduces its dependence on nif, though introduces 
explicit uncertainty due to Ai^^^. At the same time the theoretical moment is rather 
sensitive to the value of 

Let us discuss the choice of the moment in expression (^) for the 6-quark pole mass 
in more detail. It is well known that for the high moments the ordinary perturbative 
expansion in the strong coupling constant is not applicable and the Coulomb resum- 
mation is required. Recently the expansion of the theoretical moments around the 
exact Coulomb solution has been obtained in the NNL approximation p, |1^, |lT 



In our analysis we use the analytical expression of the polarization function near 
threshold obtained in refs. ||, || . The explicit formulae are too large to be presented 
here. 

The experimental moment entering eq. @ is given by the integral of the spec- 
tral density eq. (^) and is mostly saturated with the contribution of the first six T 
resonances that leads to the formula 

The leptonic widths F^, and masses {k = 1 ... 6) of the resonances are known with 



good accuracy |T2[, the electromagnetic coupling constant is renormalized to the 
energy of order mj, with the result aQEoi'^^b) = l.OTo;^ [|12|. The rest of the spectrum 
beyond the resonance region for energies larger than sq ~ (11.2 GeV)^ (continuum 
contribution) is approximated by the theoretical spectral density multiplied by the 
parameter 0.5 < t < 1.5 which accounts for the uncertainty in the experimental data 
in this energy region. Note that we assume rather large uncertainty of the continuum 
to be on the safe side however its contribution is essentially suppressed in comparison 



with the resonance one and the resulting error of the whole quantity in eq. is of 
the same order as the uncertainties introduced by the resonance contribution. 

Note that mf (or mb) can be formally extracted from an arbitrary moment of 
the spectral density. The region of allowed n however is quite restricted. Indeed, 
the low moments cannot be used in sum rules because the experimental spectrum 
is well known experimentally only for energies close to threshold due to existence of 
sharp resonances while the contribution of the continuum to these low moments is 
large in comparison with the resonance contribution. On the other hand for n > 10 
the perturbative expansion of the moments around the Coulomb solution is strongly 
divergent. This can be considered as a signal that for large n the Coulomb solution 



is not the best zero order approximation [jT3[. So n = 5 seems to be the natural and 



optimal choice. In any case the result of calculation is almost insensitive to the local 
variation of n around this value. 

From eqs. (|l],H) we find for the mixing angle \Vcb\ 



where the functions 



K{as,mb) = . = (7) 

Cr(«.)v/Mf («„mfe) 

and F{ml/ml) accumulate theoretical information depending on m^, and a^. 

Function F{ml/ml) introduces rather large theoretical uncertainty if masses of b- 

and c-quarks are considered as independent variables. However there is almost model 

independent constraint of the form 

(mfe - mj(l + 0(l/mfe,J) = fhB-mD = 3.34 GeV (8) 

where = 5.31 GeV, rh£, = 1.97 GeV denote the spin-average meson masses, 
e.g. ffiB = ^{rriB + SniB*)- We take the value of the nonperturbative corrections to 
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eq. (^) given in ref. [|T^ and use mi — mc = 3.47 GeV as the central value for numerical 
estimates. With this constraint the function F{ml/ml) becomes a function of a single 
variable F{mb). Note that in such a setting the rrib dependence of the function F{mi)) 
partly cancels the large ml dependence of the width. Furthermore using the relation 
(0) to express the 6-quark pole mass in the argument of the function F{mb) in terms 
of the moments of the spectral density results in only logarithmic dependence of the 
right hand side of eq. (|^) on mj,. 

Now one can analyze the theoretical factor K{as,mb)/F {mj,) numerically order 
by order in a^. The result reads 

f K{as,mb) y ^ 1.234(1 + 0.102 + 0.014). (9) 

V F'inib) J 

In the numerical analysis of extracting |V"c6| we use the "world average" value of the 
strong coupling constant as{Mz) = 0.118 |jl2[ and the normalization point fi ^ nib 
in the expression for the theoretical moment^. The value of 6-quark pole mass in the 
fixed order in as with Coulomb resummation is found from eq. (0). For comparison, 
the perturbative series for nib that follows from eq. (@) and the series for Cr are 

mb = 4.75(1 - 0.014 + 0.022), 

Cr = 1 -0.146 -0.064. 

Thus we find that in the expansions of the theoretical moment (as well as nib itself) 
and width expressed in terms of rrib the NNL corrections are of the order of the next- 



^ For the numerical estimates it is important to fix the aUowed range for the normahzation point 
which is present in the exphcit formula of the polarization function. The naive estimate of the 
"physical scale" of the problem /i ~ mbas ||, |l^ is not acceptable since the direct calculation of 
the NNL corrections shows that the perturbation theory for the moments blows up there |||, ^ . The 
relative weight of the NNL order corrections is stabilized at ^ ^ mi, jl^ so in our opinion there is 
no reason to use the lower normalization scales. 
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to-leading onesQ while the perturbative series for the mixing angle, or the theoretical 
coefficient eq. converges much better. 

As for numerics, we use the following central values for our experimental inputs 
(see 



121 for more detail): 

BR{B XM) = 10.5%, tb = 1.55 ps. 



Ms'"^ = 4.51 X 10-^ GeV-^ 
With these numbers we obtain the value of the matrix element \Vcb\ 

,K.,^o.o«3(52<£^)^(if^)% 



^ Mr - 4.51 X 10-4 GeV-n / aJM^ - 0.118\ , , , 

1 + 0.02^^ — E 1 - 0.01 ^ ^ 1 ± Anr,t] 

0.20 X 10-4 GeV-^ I 0.006 ' ^ 



where the nonperturbative corrections are included according to ref. |Q and A„pt ~ 
0.02 is the uncertainty in the nonperturbative contribution coming mainly from the 
uncertainty in the HQET relation between and rric eq. (^. The typical scale of 
uncertainty of key parameters is also indicated. Another important source of the 
uncertainty is the scale dependence of the theoretical moment and the experimental 
errors in the value of as{Mz) because of rather high sensitivity of the theoretical mo- 
ment to ag. In fact these uncertainties are closely related since the scale dependence 
of 7W|'^(as(/i), /i) is mainly due to the scale dependence of agifJ^) while the explicit 
dependence on fi is rather weak. The pointed error bars roughly correspond to the 
interval fi = rrib ± 1 GeV at fixed as{Mz) = 0.118. Note that the result is almost 
insensitive to the specific number of the moment used for the estimate. For example, 
the value of \Vcb\ matrix element changes approximately by 0.2% when the moment 
number changes in the interval 5 < ri < 15. 



^ There is a hypothesis that this fact is a consequence of the asymptotic character of the series 
which leads to the intrinsic ambiguity in the heavy quark pole mass fl^ . 
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Now the advantage of our approach becomes clear - large errors due to uncertainty 
in rrih is now partly shifted to more direct experimental data. Furthermore the terms 
in the perturbative expansion @ decrease rapidly which is a solid indication that the 
higher orders corrections to the obtained result are small enough. 

The main part of the experimental uncertainty is related to the uncertainty in 
the experimentally measured inclusive semileptonic width. The uncertainty in M|^^ 
comes mainly from the continuum contribution to the moments above Sq and the 
uncertainties in F^. Experimental situation is rather dynamic and data are improving 



fast that means that the experimental uncertainties will be smaller (see e.g [15|). 

Our result for the central value of the parameter |V"cf,| of the CKM mixing matrix 
eq. ( p!OD is in a good agreement with the previous estimate \Vcb\ = 0.0419 0. This 
value, however, is somewhat larger than the estimate \Vcb\ = 0.039 ± 0.002 of ref. [0. 
There is no much hope to reduce the uncertainty in the nonperturbative contribution. 
Thus the model independent result presented in the paper provides one with the most 
reliable and accurate estimate of the CKM matrix element \Vcb\ from the inclusive 
5-meson semileptonic width. 
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